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Abstract 

We briefly discuss the collinear factorization formula for the associated production of one particle and a Drell-Yan 
pair in hadronic collisions. We outline possible applications of the results to three different research areas. 
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1. Introduction 

The description of particle production in hadronic col- 
lisions is interesting and challenging in many aspects. 
Perturbation theory can be applied whenever a suffi- 
ciently hard scale characterizes the scattering process. 
The comparison of early LHC charged particle spec- 
tra with next-to-leading order perturbative QCD predic- 
tions [1] shows that the theory offers a rather good de- 
scription of data at sufficiently high hadronic transverse 
momentum, of the order of a few GeV. For inelastic 
scattering processes at even lower transverse momen- 
tum, the theoretical description in terms of perturba- 
tive QCD breaks down since partonic matrix elements 
diverge as the transverse momenta of final state par- 
ton vanish. The situation is further complicated by the 
occurence of additional parton-parton interactions yfl 
accompanying the main hard process. This effects is 
possibly ascribed to the composite and extended struc- 
ture of colliding hadrons. The nature of the latter phe- 
nomena, which at present require some kind of mod- 
elization, could be more precisely investigated if an im- 
proved theoretical description of the main hard process 
were provided. A well suited process for these kind of 
studies is represented by the associated production of 
one particle and a Drell-Yan pair, H1+H2 — > H + y* +X. 
While the high mass lepton pair constitutes the perturba- 
tive trigger which guarantees the applicability of pertur- 
bative QCD, the detected hadron H could then be used, 
without any phase space restriction, as a local probe to 



investigate particle production mechanisms. Within this 
particular process we will review how to improve the 
theoretical description to make it capable to handle the 
aforementioned collinear singularities associated to fi- 
nal state parton production at vanishing transverse mo- 
mentum. We further note that this process is the sin- 
gle particle counterpart of electroweak -boson plus jets 
associated production [3], presently calculated at nex- 
to-leading order accuracy with up to three jets in the 
final state [4]. One virtue of jet requirement is that it in- 
deed avoids the introduction of fragmentation functions 
to model the final state, which are instead one of the 
basic ingredients entering our formalism. At variance 
with our case, however, jet reconstructions at very low 
transverse momentum starts to be challenging ||5|] and 
it makes difficult the study of this interesting portion of 
the produced particle spectrum. 
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2. Collinear facrotization formula 

The associated production of a particle and a Drell- 
Yan pair in term of partonic degrees of freedom starts 
at 0(a s ). One of the contributing diagrams is de- 
picted in Fig. ([TJ. Assuming that the hadronic cross- 
sections admit a factorization in term of long distance 
non-perturbative distributions and short distance per- 
turbative calculable matrix elements, predictions based 
on perturbative QCD are obtained convoluting the rele- 
vant partonic sub-process cross-sections, da^., with par- 
ton distribution functions, /* and fj, and fragmentation 
functions, D H ^ k . The hadronic cross-sections, at center 
of mass energy squared S = (Hi +H2) 2 , can be symbol- 
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Figure 1: Example of diagram contributing to hadron production in 
the central fragmentation region to order 0(a s ) in eq.(TJ. 



Figure 2: Parton model formula, eq.|2}, for the associated production 
of a particle and a Drell-Yan pair. 



ically written as 
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where the convolution are over the momentum fractions 
of the incoming and outgoing partons and z is the en- 
ergy of the observed hadron H scaled down by the beam 
energy, yS/2, in the hadronic centre of mass system. 
The invariant mass of the virtual photon is indicated by 
Q 2 . The partonic indeces i, j and k in the sum run on 
the available partonic sub-process. The superscripts la- 
bel the incoming hadrons and the presence of crossed 
term is understood. Within this production mechanism, 
the observed hadron H is generated by the fragmen- 
tation of the final state parton k, and for this reason 
we address it as central. The amplitudes squared 01, 
however, are singular when the transverse momentum 
of the final state parton vanishes. In such configura- 
tions, the parent parton k of the observed hadron H is 
collinear the incoming parton p\ or p2- Perturbation 
theory looses its predictivity as these phase space region 
are approached. The same pattern of collinear singular- 
ities are found also in an analogue calculation in Deep 
Inelastic Scattering [8]. In both processes such singular- 
ities can not be treated with the usual renormalization 
procedure which amounts to reabsorb collinear diver- 
gences into a redefinition of bare parton and fragmen- 
tation functions. Such singularities will appear in ev- 
ery fixed order calculation in the same kinematical lim- 
its spoiling the convergence of the perturbative series. 
One possible solution is to regularize the partonic cross- 
sections with a lower transverse momentum cut-off J^l . 
It would be, however, highly desiderable to develop a 
technique to resum such logarithmic contributions to all 
orders in perturbation theory. As shown in Ref. fl, a 
proper treatment of such singularities requires the intro- 
duction of new non-perturbative distributions, fracture 
functions ifioll . M H ^ Hi (x,z), which parametrize hadron 
production at vanishing transverse momentum. They 



express the conditional probability to find a parton i 
entering the hard scattering while an hadron H is pro- 
duced with fractional momentum z in the target frag- 
mentation region of the incoming hadron. By using the 
same generalized collinear subtraction procedure devel- 
oped for fracture functions in Ref. 18J, it is then possi- 
ble to show that additional singularities discussed above 
can be properly factorized into the redefinition of bare 
fracture functions |g]. The fact that the collinear sub- 
traction used in Drell-Yan process is the same as in 
DIS supports both the universality of collinear singu- 
larities among these processes with a different number 
of hadrons in the initial state and the factorization for- 
mula in eq. Fracture functions inhomogeneous evo- 
lution equations 1 10], which can be readily derived start- 
ing from this generalized collinear subtraction, can then 
be used to resum to all orders such type of collinear 
singularities. The use of fracture functions opens the 
possibility to have particle production already to 0(a° s ), 
since the hadron H can be non-pertubatively produced 
by a fracture function M itself. The lowest order parton 
model formula can be symbolically written as 
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where the superscripts indicate from which incoming 
hadron, Hi or H2, the outgoing hadron H is produced 
through a fracture functions. This production mecha- 
nism is sketched in Fig. (0. So far we have only consid- 
ered 0(a s ) corrections in the central fragmentation re- 
gion, eq. ([TJ, to the parton model formula, eq. (fj) . In or- 
der to complete the calculation to 0(a s ) we should also 
consider higher order corrections in process initiated by 
a fracture functions. In this case, in fact, the hadron H 
is already produced by these distributions and therefore 
final state parton in real emission diagrams should be 
integrated over and results added to virtual corrections. 
One of the contributing diagrams is depicted in Fig. O). 
The general structure of these terms is 
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Figure 3: Example of diagram contributing to 0(a s ) corrections in the 
target fragmentation region, eq.(3J- 



d(r H,T,(l) 



dQ 2 dz 2n 

We refer to this corrections term as to the target frag- 
mentation contribution. The calculation is, a part from 
different kinematics, completely analogue to the inclu- 
sive Drell-Yan case. While soft singularities cancel 
among real and virtual diagrams, the collinear singu- 
larities appearing in this term are renormalized by the 
homogenous renormalization term in bare fracture func- 
tions and by bare parton distributions functions. Adding 
all the various contributions, the resulting p, -integrated 
cross-sections is then infrared finite [6] and can be sin- 
tetically written as 



dcr H 
dQ 2 dz 
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We refer to the previous equation as to the collinear 
factorization formula for the process under study. The 
next-to-leading coefficients C' ; and K'J have been cal- 
culated 111 111 , making the whole calculation ready for 
numerical implementation. In the following section we 
will try to outline some applications of the proposed for- 
malism. 



3. Applications 

a) The first possible benchmark process for the calcula- 
tion could be strangeness production associated with a 
Drell-Yan pair, p + p — > V + y* + X, where V gener- 
ically indicates a A or A hyperon. The first process, 
according to so called leading particle effect, should be 
sensitive to the spectator system fragmentation into A 
hyperon at very low transverse momentum and there- 
fore to be modelled with the help of fracture functions 
in dcr H ' T . In the anti-hyperon production case no such 
effect should be present so that A are expected to be 



mainly produced by the fragmentation of final state par- 
ton as described by dcr H ' c . The measurament of such 
process for different Q 2 could allow to test the frac- 
ture functions scale dependence embodied in their pe- 
culiar evolutions equations and the validity of the pro- 
posed factorization formula in all its components. With 
this respect, a comparison of strange particle produc- 
tion in hadronic collisions and DIS would be extremely 
interesting since longitudinal momentum spectra of A 
and A hyperon has been measured quite accurately in 
charged current neutrino DIS Hill , both in the current 
as well in the target fragmentation region. 

b) The formalism may found application in the study 
of single diffractive hard process, p + p — > p + y* + X, 
where the outgoing proton has almost the incoming pro- 
ton energy and extremely low transverse momentum 
with respect to the collision axis. This process has been 
intensively analyzed in the DIS at HERA, revealing its 
leading twist nature. From scaling violations of the 
diffractive structure functions 111 311 and dijet production 
in the final state Ill4ill5ll quite precise diffractive parton 
distributions functions have been extracted from data, 
which parametrize the parton content of the color singlet 
exchanged in the f-channel. The predictions for single 
diffrative hard processes measured at Tevatron, based on 
diffractive parton distributions measured at HERA and 
assuming factorization, has indeed revealed that these 
processes are significantly suppressed in hadronic colli- 
sions, see the very recent analysis reported in Ref. [16]. 
Recalling that these distributions are fracture functions 
in the z — > 1 limit, the present formalism can then ap- 
plied to next-to-leading order accuracy, the main contri- 
bution to the cross-sections coming from the target term, 
dcr H ' T . Such term can be eventually recast in triple dif- 
ferential form in xjp ^ 1 - z, virtual photon rapidity, y, 
and invariant mass Q 2 and evaluated at next-to-leading 
order by using the appropriate coefficient functions 01711 . 
In this way factorization tests could be porformed at 
fixed xjp to avoid any Regge factorization assumption 
on DPDF while the y dependence, giving direct access 
to the fractional parton momentum in the diffractive ex- 
change, p, allows to test factorization in a kinematic re- 
gion which avoids DPDF extrapolation. Finally, the Q 2 
dependence of the cross-sections could be used to in- 
vestigate how factorization breaking effects eventually 
evolve with the hardness of the probe and to which ex- 
tent the factorized formula M ® / actually works. 

c) The calculation has been performed to make 
predictions for cross-sections integrated over partonic 
transverse momentum. To this end, divergent contribu- 
tions, related to parton emissions at vaninshing trans- 
verse momentum, are factorized into fracture functions. 



3 



In the case, however, that cross-sections are measured 
down to a minimum but still perturbative hadronic trans- 
verse momentum, the latter constitutes a natural in- 
frared regulator for the partonic matrix elements. The 
central production term, do JI,c,< - 1 \ can be used to esti- 
mate hadron production as the fragmentation process, 
parametrized by fragmentation functions, were happen- 
ing in the QCD vacuum. The hadronic cross-sections 
can be recast in a triple differential form in Q 2 , produced 
hadron transverse momentum p, and pseudo-rapidity 77 
to predict charged particle spectra or multiplicity. A par- 
ticular interesting observable which can also be recon- 
structed is the differential cross-sections differential in 
cos </>, where <p is the angle formed by the virtual pho- 
ton and the detected hadron in the center of mass sys- 
tem. This observable has been shown to be sensitive to 
the contamination of the so-called underlying event 111 80 
to jet observable and has been used also to investigate 
underlying event properties in Drell-Yan process II 1 911 . 
Although predictions made on the present formalism 
take into account the radiation accompanying one single 
hard scattering per proton-proton interactions, it can be 
nevertheless used as a reference cross-sections to gauge 
the impact of new phenomena. The dcr H ' c term, altough 
formally 0(a s ), is a tree level predictions and possible 
large higher order corrections may be expected espe- 
cially in the forward region at large transverse momen- 
tum, as already found in DIS analysis l20ll . We finally 
note that the recent analysis of underlying event struc- 
ture in Drell-Yan events is performed in a wide elec- 
troweak boson mass window [19]. We emphasize in- 
stead that more insight on the underlying event dynam- 
ics could be accessed if cross-sections could be mea- 
sured, as in the single diffractive case discussed above, 
at different virtual photon virtualities Q 2 and transverse 
momentum qj. 



4. Conclusions 

We have briefly reviewed a perturbative approach to sin- 
gle particle production associated with a Drell-Yan pair 
in hadronic collisions. On the theoretical side we have 
shown that the introduction of new non-perturbative dis- 
tributions, fracture functions, allows a consistent factor- 
ization of new class of collinear singularities stemming 
for configurations in which the parent parton of the ob- 
served hadron is collinear to the incoming parton. The 
scale dependence induced by this generalized factoriza- 
tion is driven by Altarelli-Parisi inhomogeneous evolu- 
tion equations for fracture functions which allow the re- 
summation to all orders of this new class of collinear 
logarithms. The factorization procedure does coincide 



with the one used in DIS confirming, as expected, the 
universal structure of collinear singularities among dif- 
ferent hadron initiated processes and supporting the 
collinear factorization formula proposed in eq. On 
the phenomenological side we have briefly discussed a 
few applications in which different aspects of the for- 
malism could be tested and compared towards results 
already obtained in the DIS target fragmentation region. 
A good theoretical control on the perturbative compo- 
nent indeed allows the investigation on new phenomena 
appearing in hadronic collisions, for example the rapid- 
ity gap probability suppression in hard diffractive pro- 
cesses with respect to diffractive DIS and the investiga- 
tion, although indirect, of multiple parton-parton inter- 
actions. 
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